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Abstract. An analytic proof is proposed of Wiener's theorem on 
factorization of positive definite matrix-functions. 



1. Introduction. The series of papers of Wiener [15], [16], [18], and 
Helson and Lowdenslager [4], [5] led to the following 

Matrix Spectral Factorization Theorem: Let 

\z\ = 1, be a positive definite rxr matrix-function with integrable entries, 
Sjk{z) G Li(T), defined on the unit circle. If 

(2) logdet 5(2) G Li(T), 
then S{z) admits a factorization 

(3) Siz)=xHz)ixHz)r, 

where x^{z), \z\ < 1, is an analytic rxr matrix-function with entries 
from the Hardy space H2, and detx^{z) is an outer function. 

The equation (3) is assumed to hold a.e. on the unit circle T and 
{x^y = ix'^)'^ is the Hermitian conjugate of x~^- 

The spectral factorization (3) is unique up to a constant unitary mul- 
tiplier, and the unique spectral factor with additional requirement that 
X^(0) is positive definite is called canonical. 

It is well-known that the condition (2) is also necessary for the existence 
of the factorization (3). 

Spectral factorization problem was originally formulated in the scalar 
case by Wiener [14] and Kolmogorov [9] , independently from each other, 
in connection with developed linear prediction theory of stationary sto- 
chastic processes. Wiener [16] used the matrix generalization of the prob- 
lem in the study of multidimensional stationary time series. Since then 
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the spectral factorization has became an important tool in solution of 
various applied problems in Control Engineering and Communications. 

The prediction theory of stationary processes itself, namely the idea 
derived from the theory of linear least-squares method, is helpful to prove 
the spectral factorization theorem in the scalar as well in the matrix case 
and Wiener successfully exploited this idea. Helson and Lowdenslager 
derived the theorem from the theory of invariant subspaces and this proof, 
together with its generalizations to positive operator valued functions, is 
well presented in [3]. 

However, in the scalar case, the theory of the Hardy spaces Hp provides 
a simple solution to the spectral factorization problem. Furthermore, the 
spectral factor can be written in an explicit form (see [10]) 



No analog of the formula (4) exists in the matrix case since the equation 
exp{A + B) = expAexpB breaks down for non- commutative matri- 
ces and, in general, the factorization problem becomes more demanding. 
Nevertheless, one can observe that the formulation of the spectral factor- 
ization theorem in the matrix case does not go beyond the theory of Hp 
spaces and naturally its proof had been searched within this theory. In 
the present paper we provide such proof on 50 years anniversary of this 
important theorem. 

Together with the existence theorem of the spectral factorization, the 
challenging problem became the actual approximate computation of the 
spectral factor x~^{^) fo^' ^ given matrix-function (1) since, as it was 
mentioned above, this procedure found its immediate applications in so- 
lutions of various practical problems. Wiener [17] made first efforts to 
create a sound computational algorithm of matrix spectral factorization 
and since then several dozens of papers have addressed to the develop- 
ment of such algorithm. It is surely hopeless to give here a short but still 
comprehensive account on all these developments (see [7], [13]). Wiener 
and his followers [17], [11] as well as Kolmogorov's students in Russia 
[12], [20] were using methods of Functional Analysis for solution of the 
problem. A Newton-Raphson type iterative algorithm was constructed 
by G. Wilson [19]. Kalman's state-space approach to the Wiener filter- 
ing problem [8] became fruitful for spectral factorization as well, and by 
this way the spectral factorization problem was actually reduced to the 
solution of algebraic Riccati equation. Until recently this method was 
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considered as most effective for matrix spectral factorization (see [7], p. 
206), though it can be used only for rational matrix functions. Numerous 
algorithmic type improvements were proposed by the different authors 
in this direction. 

An absolutely new approach to the spectral factorization problem in 
the matrix case was developed in [6] (see also [1]). Without impos- 
ing on matrix- function (1) any additional restriction, apart from the 
necessary and sufficient condition (2), effective method of approximate 
computation of x"*" is proposed. This is the first time that the meth- 
ods of Complex Analysis and Hardy Spaces were used for construc- 
tion of the algorithm, which naturally turned out to be very efficient 
as corresponding software implementation confirms (see demo version at 
www.rmi.acnet.ge/SpFact). The decisive role of unitary matrix-functions 
in the factorization process is revealed which, by flexible manipulations, 
completely absorbs all the technical difficulties of the problem leaving 
very few and simple procedures for computation. This method of matrix 
spectral factorization is undoubtedly as simple as possible and in the 
present paper we would like to demonstrate that in fact it contains an 
analytic proof of the existence theorem as well, which was formulated in 
the beginning. No other assumptions beyond the theory of Hardy spaces 
are used for this purpose. 



2. Notation. Lp{T), p > 0, is the class of p-integrable complex func- 
tions defined on the unit circle, and Hp denotes the Hardy space of ana- 
lytic functions in the unit disk, 

Hp^lfe A{D) : sup / " |/(re^*)r < oo] . 

I r<l Jo J 

A function Q{z) e Hp is called outer, we denote Q{z) e H^, if 

Q{z) = c . exp Q- ^ log |Q(e^*)| d?j , 
where |c| = 1. 

The nth Fourier coefficient of / G ^i(T) is denoted by c„(/) and, for 
p > 1, L+(T) := {/ e Lp{T) : c„(/) = forn < 0}. The spaces L+(T) 
and Hp are naturally identified, so that we can speak about the value of 
function / e L+(T) in z e D. 
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Let Lp-{T), N > 0, p > 1, he the set of functions / from Lp{T) for 
which c„(/) = whenever n < —N, and let Lj^j be the set of analytic 
polynomials whose nonzero coefficients range from to N. 

If M is a matrix (resp. matrix-function), then M denotes the matrix 

(resp. matrix-function) with conjugate entries and M* :— . The 
upper-left mxm submatrix of M is denoted by Ml"*'"*' . 

We say that matrix-functions have some property, say, belong to Lp{T) 
or are convergent, etc, if their entries have this property. 

A rx r matrix U is called unitary if UU* = Ir, where is the identity 
matrix of dimension r. Since the rows and columns of U are orthonormal, 
its entries are bounded by 1. A unitary matrix-function U{z) means that 
it is unitary for a.a. z 

The notation diag(Mi,M2, • • • • "r) stands for the diagonal rxr matrix 
with corresponding entries on the main diagonal. 

3. The uniqueness of spectral factorization. In the proof of the 
convergence property of the above mentioned algorithm given in [6], at 
least formally, the existence of spectral factorization is used. However, 
one can observe that rather the uniqueness than the existence of spectral 
factorization provides this convergence. So we start with a simple proof 
of the uniqueness theorem (cf. [2]), emphasizing that it can be obtained 
without a priori knowledge of the existence of the spectral factorization 
itself. 

We use the following generalization of Smirnov's theorem concern- 
ing functions from the Hardy spaces Hp (see [10], p. 109): Let f{z) — 
g{z)/h{z), where g e Hp^, pi > 0, and h e Hp^, p2 > 0. // the boundary 
values /(e**) e Lp{T), p > 0, then f E Hp . 

Uniqueness Theorem: // 

(5) S{z)^xti^)ixti^))\ i = l,2, 

are two spectral factorizations of a given spectral density S{z), then 

(6) xti^) = xti^)-U, \z\<l, 
for some constant unitary matrix U. 

Proof. It follows from (5) that Xi{^){Xi{^))* = X2 (^)(X2 (-2))* a.e. 
on T, so that 

(7) {Xt{z))-'xt{z){{xt{z))-'xtiz)y = for a.a. z e T. 
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Thus the analytic matrix-function 

(8) U{z):={xt{z)r'xti^), kl<l, 

is unitary on the boundary for e T. Consequently U{e'^) G L^{T). 

Since j = 1,2, are spectral factors, it is assumed that their 

entries are from H2 and detx^(2;), j — 1, 2, are outer analytic functions. 



so that entries of 



U{z) 



1 



det xtiz) 

can be represented as ratios of two functions from i/2/r and H^,^,, re- 
spectively. Hence, we can use the generalization of Smirnov's theorem to 
conclude that U{z) e H^, i.e. C/(e**) G L^{T). 

By changing the roles of Xi and xt this discussion, we get 

ixtizT'xtiz) e H^. 

But iU{z))* = {Xiiz))~'^X2i^) for a.a. ^ G T, by virtue of (7). Thus, we 
have 

f/(e^*) e L+ (T) andf7(^eL+(T), 

which implies that (8) is a constant matrix-function and (6) follows. □ 
4. Main Lemmas. 

Lemma 1. For anymxm matrix-function Fm{z) of the form 

•••0 \ 



(9) Frr^iz) 



( 1 








1 








1 













1 



0- 

\pi{z) (p2iz) (fsiz) ■ 

\z\ = 1, where f^ G H!^ and (pj G (T), j = 1,2, . . . ,m — 1, for some 

positive integer N, there exists a unitary matrix-function Um{z) of the 
form 

( uix{z) uviiz) ■■■ uimiz) \ 

U2l{z) U22{z) ■■■ U2m{z) 



(10) 



Um{z) 



^-^m— 1,1 (^) 



^m—l,m (^) 



\ Uml{z) Um2{z) 



Ur 
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\z\ = 1, where 

Ujk G -f'[jv] ) i) ^ = 1) 2, . . . , m, 

and 

(11) deiUm{z) = l, 
\z\ — 1, such that 

(12) F^{z)Um{z)eL+{T). 

For two dimensional matrices this lemma is proved in [6], and it is 
generalized for any dimensional matrices in [1]. We emphasize that the 
lemma can be proved without any reference to the existence theorem of 
matrix spectral factorization. Furthermore, a system of linear equations 
which provides the coefficients of functions Ujk-, j, k = 1,2, . . . ,m, when- 
ever N negative coefficients of (pj, j = 1, 2, . . . , m — 1, and N positive 
coefficients of /"*" are given, can be written and solved explicitly (see [1], 
p. 22). 

Lemma 2. For any mxm matrix-function Fm{z), \z\ = 1, of the form 
(9), where /"*" e H2 and cpj e L2{T), j = 1,2, . . . ,m — 1, there exists a 
unitary matrix-function Ujn{z) of the form (10) satisfying (11) a.e. such 
that 

(13) ujkiz) e L^iT), j,k = 1,2,..., m, 
and (12) holds. 

Proof. Let f!^'' (z) be the L2-approximation of matrix-function Fm{z) 
where the entries (fij{z), j — 1, 2, . . . , m — 1, are approximated by there 
Fourier series 

00 

n=-N 

Then we can use Lemma 1 which provides the existence of unitary matrix- 
function Um\z), 



(14) 



UL''\z){U^''\z)r ^ Im (a.e.). 
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such that 
(15) 





fu^^iz) u[?{z) ... 


u[^iz) \ 






4m (^) 




'^m-l,l(^) '^m-l,2(^) ■ ■ ■ 


u^^^ (z) 






Ummi^Z^ J 


(16) 


det[/(f)(z) = l 


(a.e.), 


and 






(17) 


FL''\z)Ul^\z)eLt{T). 



7/(^) P r+ 



Now, a convergent subsequence can be extracted from {uln^}'^^^. Fur- 
thermore, if we require in addition, say, positive definiteness of the matrix 

F^n'^U'^n \{]) which can be achieved by multiplying, if necessary, u!h^\z) 
from the right by a constant unitary matrix with determinant 1, then 
u!h^^ {z) itself converges at least in measure, 

(18) U^\z)^UUz), 

as N ^ oo. These facts were proved in [6] for two dimensional case and, 
in the similar way, this can be done for any dimensional matrices as soon 
as the explicit form of Um \z) is obtained. Anyway, the convergence 

(18) , together with boundedness of unitary matrix- functions, guaranties 
that we can pass to the limit in (14)- (17), so that unitary matrix-function 
Um{z) in (18) satisfies the desired conditions (10)-(13). □ 

5. Proof of the theorem. First perform the lower-upper triangular 
factorization of (1) with positive entries on the diagonal, i.e. take 

(19) S{z) = A{z){A{z)y for a.a. ^ e T, 
where 

... 

f22{z) ■■■ 



A{z) 



//ii(^) 

f2l{z) 



\ 



and /- 



33 \ 



\frl{z) 

> for a.a. z G T, j 



fr2(z) ■ ■ ■ frr(z)J 

= 1,2, ... ,r. This factorization can be 



achieved pointwise by the well known theorem of Linear Algebra. 
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Since Sjj{z) G Li{T), j = 1, 2, . . . , r, by virtue of equations 

k=l 

(see (19)), all the entries of A{z) are square integrable, 

fjk{z) e L2(T), 
^ j f, ^ < j. Furthermore, since 

r r ^ 

5^1og^(^) = logJ]^(^) = logdetyl(;2) = -logdet5(;2) G Li(T) 

3=1 j=l 

(see (2)), we have 

(20) log/,,(2;)GLi(T), j = l,2,...,r. 

The condition (20) provides that 



= exp (^^ 1^ log /,,(e^*) j , J = 1, 2, ... , r. 

arc outer analytic functions satisfying \ f^{z)\ = fjj{z) for db.db. z G T. 

If wc denote the ratio fj^{z)/fjj{z) by Uj{z), then |Mj(-2)| = 1 for a. a. 
z & T, j — 1,2, . . . ,r, and U{z) — dieig{ui{z) , U2{z) , . . . , Ur{z)) is unitary 
matrix-function. Thus, for matrix-function M{z) — A{z)U{z), we have 

(21) S(z) = M(z){M(z))* for a.a. zeT, 

where 

/ fH^) ••• \ 

<^2l(^) /2^(2) 



(22) M(^) 











ipr-l,l{z) ipr-l,2{^) 



frUiz) 
(Pr,r-l{z) fr{z)) 



f+{z) G HO, G L2(T), 1 < J < r, 1 < A; < J. Define 

(23) Miiz) = M{z) said M^{z) = Mm-i{z)VUz), m = 2,3,...,r, 
where K?i(^), m = 2, 3, . . . , r, are unitary matrix-functions 

(24) Vm{z)V:^{z) = Im, for a.a. z G T, 

constructed recurrently as follows: let Fm{z) be the mxm matrix of the 
form (9), where its last row coincides with the last row of Mj^Z^\z), and 
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let Um{z) be the corresponding unitary matrix-function which existence 
is proved in Lemma 2, so that (11) and (12) hold. Define Vm{z) as the 
block matrix-function 



(25) Vm{z) 



Um{z) 
-^i — m 



m = 2,3, . . . ,r. It is assumed that Vr{z) — Ur{z). Obviously Vm{z) is 
unitary matrix-function and (sec (11)) 

(26) det Vm{z) = det Um{z) = 1 for a.a. z eT. 
Pay attention that the following equation holds 

(27) {M^::-!^"^-"^ (^) 0^ . F^{z) = M^-f {z) , 
while 

(28) M'^::Z,''^-'\z)eLt{T) 

lor m — 2 (see (23), (22)) and for m > 2 as well, because of the con- 
struction process (see (30) below). Indeed, by virtue of (23) and (25), we 
have 

(29) M|r-](^)=Mi"!:7l(^)C/m(^), m = 2,3,...,r. 
Thus, it follows from (29), (27), (28), and (12) that 

(30) Mi-'-l {z) = (^^^Ti'^""'^ (^) 0^ Frr.{z)Um{z) e L+(T). 
For m = r, we have 

(31) Mr{z) e L+{T), 

and we conclude that Mj.{z) is a spectral factor of >S'(^), 

(32) x-'iz) = Mr{z). 
Indeed, since (21), (23), and (24) hold, we have 

S(z) = Mr(z)(Mr(z))* for a.a. zeT, 

and it remains to show that det Mr{z) is outer, where Mr{z) is assumed 
extended in the unit disk D in this case, by virtue of (31). 
The equations in (23) imply that 

(33) Mr{z) = M{z)U2{z)U3{z) ■ ...Ur{z) for a.a. ^ e T. 
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Hence, taking into account (33), (22), and (29), we have 

(34) det Mr{z) = fi{z)f^{z) . . . f^{z) for a.a. zeT, 

The both sides of (34) are functions from H2/r and they coincide on the 
boundary almost everywhere. Thus the equation in (34) is vahd inside the 
unit circle for each z E D, and since each fj~{z) is outer, j = 1, 2, . . . , r, 
their product det Mr{z) is outer as well. 
The proof of the relation (32) is completed. 
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